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We consider a wide class of increasing Levy processes perturbed by an independent Brownian motion as a degradation model. 
Such family contains almost all classical degradation models considered in the literature. Classically failure time associated 
to such model is defined as the hitting time or the first-passage time of a fixed level. Since sample paths are not in general 
increasing, we consider also the last-passage time as the failure time following a recent work by Barker and Newby Mm. We 
address here the problem of determining the distribution of the first-passage time and of the last-passage time. In the last section 
we consider a maintenance policy for such models. 
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For several decades, degradation data have been more and more used to understand ageing of a device, instead of 
only failure data. The most widely used stochastic processes for degradation models belong to the class of Levy 
processes. More precisely, the three main models are the following ones: (a) Brownian motion with (positive) 
drift; (b) gamma processes; (c) compound Poisson processes. More generally we consider a broad class of Levy 
processes corresponding to subordinators perturbed by an independent Brownian motion: 



where {Gt, t > 0} is a subordinator, i.e. a Levy process with non decreasing sample paths. Since jumps of 
{D t , t > 0} are issued from {Gt, t > 0} and are positive, we recall that we say that {D t , t > 0} is spectrally 
positive. This process can be characterized in terms of Levy exponents: 



Exponent (f>s is associated to the Brownian motion and 4>q to Gt, which is in all generality a jump process. 
It follows that the Levy measure of {D t , t > 0} is the same as that of {Gt, t > 0} that we will denote by 
vd(<1x) = Q(dx). Furthermore we will suppose that measure Q{.) admits a density with respect to the Lebesgue 
measure, i.e. that Q(dx) — q(x)dx for some density q(.). In the following we will also need 



i.e. (pd{u) is such that E[e uDt ] = e tipD ^ . We recall, since {Gt, t > 0} is a subordinator, that may write in this 
case ipd{u) in the following way 



1 Introduction and Model 



Vi > , D t = G t + aB t 




R\{0} 



<Pd(u) = 4>d{iu) = <pg(u) + 7rii 2 <7 2 , 




\e 



—ux 



l]Q(dx) + -u 2 a 2 , 
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for some /i > 0. We consider in this paper several approaches for modelling degradation of a device and its failure 
time. Failure time can traditionally be derived from a degradation model by considering the first hitting time Tf, 
of a critical level b > 0. The first-passage time distribution has been already derived for the particular case of two 
sub-models. In the case of Brownian motion with drift (corresponding to Gt = [it, /j, > 0), it is the well-known 
inverse Gaussian distribution, see Jl5| ] for instance. For the pure gamma process (i.e. a = and {Gt, t > 0} 
is a gamma process), it has been studied by Park and Padgett [H. Moreover they proposed an approximation 
for the cumulative distribution function of the hitting time based on Birnbaum- Saunders and inverse Gaussian 
distributions. 

Recently a new approach to define the failure time was proposed by Barker and Newby [Q] that consists in 
considering the last passage time of degradation process {D t , t > 0} above b. As explained in that paper, this is 
motivated by the fact that, even if {Dt, t > 0} reaches and goes beyond b, resulting in a temporarily "degraded" 
state of the device, it can still always recover by getting back below b provided this was not the last passage time 
through b. On the other hand, if this is the last passage time then no recovery is possible afterwards and we may 
then consider it as a "real" failure time. Of course, this discussion about modelling failure time by the first or last 
passage time becomes irrelevant whenever process {D t , t > 0} has non decreasing paths (which is not the case 
e.g. of the Brownian motion) since in that case both quantities coincide. 

In this paper we then investigate these quantities for a rather wide class of so-called perturbed process. In Section 
^ we provide the Laplace transform of the first passage time T& with penalty function involving the corresponding 
under and overshoot of the process. We then confront this approach to related recent existing results on such 
passage times distributions in the general theory of Levy processes, that introduces the notion of so-called scale 
functions. The case of several sub-models is reviewed (or revisited) : in these cases the probability distribution 
function (pdf) and/or the cumulative distribution function (cdf) can be computed explicitly, or at least numerically. 
In conclusion of this section we propose an alternative degradation process that takes into account the fact that the 
process cannot be in theory negative and suggests that {D t , t > 0} be reflected at zero. In that setting we use the 
aforementioned recent results in the theory of Levy and reflected Levy processes to obtain the joint distribution 
of the first passage time jointly to the overshoot distribution. In Section || we study the case where failure time 
corresponds to the last passage time above b and derive its distribution in the non reflected and reflected case. 
Finally we consider in Section Q a maintenance policy problem inspired by [Q] and derive distribution of related 
quantities. 

To conclude this introduction, we make precise wh ere in the pres ent paper previously published results are 



reviewed and what is actually novel. Proposition |2.2| in Section |2.1| is new, but its proof is similar to the one 



corr esponding to proof of Remark 4.1 as well as Expressions (4.4) and (4.5) of Garrido and Morales 0160. Section 



2.3| recalls facts (with short proofs) previously establishe d in the literature that are useful later on. On the other 



hand and to the best of our knowledge, Theorems [3^ and |3.3| in Section |3j concerning last passage times may be 
linked to Chiu and Yin [fujl, Baurdoux [^J and recent paper Kyprianou et al. [ |2l| ] but are otherwise genuinely new. 
Similarly Section |] deals with determining reliability quantities features unheard-of results. 

2 First-passage time as failure time 

We consider here the hitting time distribution of a fixed level b > by the perturbed process {D t , t > 0}: 

T b = inf {t > ; D t > b} 

which we remind is a.s. finite. We study below the distribution of (T&, DT b -, Dr b ) by determining the following 
quantity 

<j) w (5,b)=E(e~ 6T >w(D Tb _,D Tb )) (1) 

where S > and w(., .) is an arbitrary continuous bounded function that will be referred to as penalty function. 
In the following we will drop the subscript when there is no ambiguity on w(., .) and then write <fr(6, b) instead of 
(f> w (5, b). We then determine (|l]) in the general case and then illustrate our results to sub-models, some of which 
distribution of T), has already been obtained. 

2. 1 General case 

We are interested in the case where process {Gt, t > 0} is general. To this end, we use a well known technique 
that consists in approaching the jump part process in { Gt , t > 0} by a compound Poisson processes which, as said 
in the Introduction, is similar to the one used in [|l6| (for more details see Appendix A.l in [|6|). More precisely 
this process can be pointwise approximated by a sequence of compound Poisson processes ((S(t, n)) t >o) n eN such 
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1. (S(t, n)) nS N is increasing for all t > 0, 

2. /it + lim, i _ i . 00 n) = G t for all £ > 0, 

3. for all n, (S(t, n))t>o has intensity A n and jump size with c.d.f. P n (x) with 

A„ = Q(l/n) (2) 

P " (X) = Q(l/n) ^/"> (3) 

where Q(x) :— Q([x, +oo)). Note that Q defines measure such that Q(dx) — —Q(dx). 

Note that this approach is particularly interesting when A„ = Q(l/n) — > 0(0) = Q([0,+oo)) = +oo as 
n — > oo, i.e. when process has infinitely many jumps on any interval. Intuitively {S(t, n), t > 0} is obtained 
from {G t , t > 0} by discarding all jumps that are of size less than 1/n. Since {S(t, n), t > 0} increases towards 
{A, t > 0}, we have that 

Tb \T b , n -> oo, a.s., (4) 

where T b ™ is the hitting time of level b of the truncated process {-D™, t > 0} defined by D™ = 5(t, n) + crBt for 
any i > and any n G N. We remind that T b ™ is also a.s. finite. In fact T b n may be described as a ruin time (i.e. 
the first hitting time of of a stochastic process) in the following way: 

T b ™ = inf {t > ; b-p,t- S(t, n) + aB t < 0} 

and we are interested in the Laplace transform <j) n (5) := ¥.(e~ ST * w(Dj,„ _ , D^ n )) of T fc ra with penalty function 
w(., .) for all S > 0. Let p n = /9 n (<5) be the positive solution to the following equation: 

poo a 2 

A„ y e~ pnX dP n {x) =X n + 6- —pi + p Pn (5) 
that we will call generalized Lundberg equation. We start by showing convergence of p n as n — >• oo. 

Proposition 2.1 p„ converges as n oo to the unique solution p > to the following generalized Lundberg 
equation: 

5-yP 2 = ( Pg(p) 6 = <p D (p) (6) 

Proof: Thanks to Expressions (||) and of A n and c.d.f. P n , we may rewrite (Jfy in the following way 

f°° a 2 [°° {°° a 2 

\ e-P nX Q(dx) = Q(l/n) + 6 - — p 2 n + np n <^=> / e-" nX Q(dx) = Q{dx) + 5 - — P 2 n + p Pn 

Jl/n ^ Jl/n Jl/n ^ 

„2 

2 



s - Y Pn + PPn + 



/•OO 

/ (1 - e- pnX ) Q{dx) = 0. 

Jl/n 



Thus p n is the only positive solution to equation f n (z) = 0where/ n (.z) := 5— ^-z 2 +pz+f^J n (1 — e zx ) Q(dx). 
Let us note that (/ n )neN increasingly converges pointwise towards 

2 /-oo 

f(z) = 5-—z 2 + pz + (1 - e — *) Q(di) = 8 - tp D (z), 
* Jo 

so that p n converges increasingly towards p* := sup ngN p n . Besides one can verify that f(z) = admits an 
unique solution on (0, +oo), which is solution p to Equation (^J). Thus p* is less than or equal to solution p and we 
prove that we in fact have equality p* = p which is achieved by showing that f(p*) = 0. Indeed, using inequality 
< zx for all z, x > and since f n (pn) — 0, we have 

= |/(P*) - /n(Pn)| < \f{p*) ~ f(Pn)\ + \f(Pn) ~ fn( P n)\ 

< \f(p*) - f(Pn)\ + (1 - e-P" X ) Q(dx) < \f( P *) - f(p n )\ + Pn \ xQ(dx) 



nl/n 

< \f(p*) - f(p n )\ + P / xQ(dx) since p n < p* < p. 



(7) 
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We recall that the fact that {G t , t > 0} is a subordinator (a non decreasing Levy process) implies that / °°(1 A 

x)Q(dx) < +00 (see e.g. (2) p. 72 of [^|), hence Jq 1 ^" xQ(dx) — > 0. Remembering that / is a continuous 
function, this implies that (fy tends to zero as n — > +00, hence /(/?*)= 0. □ 

The Laplace transform 0„((5) with penalty function w(., .) of T£ is given through the following which is a 
particular case of Theorem 2 of p9| ] adapted to our context: 

Theorem 2.1 Let w{. 1 .) be a bounded continuous function and define 

w„(x) = / w(x,y - x)dP n (y). 

J X 

Thenb^f 4> n (S,b) :— ¥,(e~ STb w(D 7 ^ n _, D^ n )) satisfies the renewal equation 

4> n (6, b) = 4^(5,-)* g n (8, ■ )(&) + K (5, b) (8) 
where functions <?«(•,•) and h n (-,-) are defined by 

gn (s,y) = ^ f y e -i-2^ +P „( m v- S ) f°° e - P -m*-s) dPn{x)ds (9) 



(T 



11 



a 2 
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h n (S,y) = e -l-^/- 2 +p^)\v + ^ I' e -l-W*>+Pnmv-) I e->>"W x - s ^ n (x)dxds. (10) 
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Proof: With notations of [[29|], we have g n {5, y) expressed as in (1.10) therein with b := b(S) = — 2/i/cr 2 + p n {8), 
A := A„, P(-) :— P n {-) and D = a 2 / 2. Still with notations of p9[], and in Theorem 2 therein, we see that function 
y 1 ^ /i„(<5, y) is the sum of e~I~ 2A1 / cr +' 3 ™(< 5 )]2' and some function g w (-) defined in Expression (2.8) of [ j^ that 
depends on u> n . It is easy to verify that this function is the last term on the right-handside of dl(i|). □ 

Passing on the limit n — > +00 in Theorem 2.1 yields the following renewal equation for function ([!]): 



Proposition 2.2 Let uj(x) := w(x, y — x)Q(dy). Function (f>(S, •) = (f> w (S, •) satisfies the renewal equation 

c/>(5,b) = <t>(S,-)*9(S,-)(b) + h(5,b) (11) 
where functions <?(•,•) and •) are defined by 

g (S,y) = 4 / e -[- 2 ^ 2 +"WK^) / e-fW^-'JQCdsJda (12) 



a 2 



(13) 



h(S,y) = e -[- 2 »/° 2 +p{S)\y + JL j e -[- 2 »/° 2 +p(my-s) j e -P^- s ^{x)dxds. 

Hence (f)(S 7 b) is given by the Pollaczek-Kinchine like formula 

00 

<t>(S,b) = J29* k (S,.)*h(6,.)(6,b). (14) 

fe=0 

Note that ( |l4| ) is analogous to Expression (4.2) in [n6p. 

Proof: Let us prove that A„w„ converges to lu. This is easily seen by remembering that A„ = Q(l/n) and thus 
that, by (ph, 



X n u n (x) = - j w(x, y - x)l {y > 1/n }dQ(y) 

which converges to the desired expression, remembering that —dQ(y) = dQ(y). Convergence of h n to h follows 
from (|To|). In the same way, \ n e~ Pn ^ s ^ x ~ s ^dP n (x) converges to e~ p ^ s ^ x ~ s ^Q(dx), yielding convergence 
of g n to g thanks to (||). □ 

As announced in the Introduction, it is also possible to use the theory of Levy processes to propose a differ- 
ent approach for determining the joint distribution of the hitting time Xj, jointly to the state of Dq\, using scale 



rfTrih. 
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Proposition 2.3 (Theorem 1 (4) [fo|]) Let us define for all 5 > the scale function W^ s \ through its Laplace 
transform, and by 

e - Xx W (5) (x)dx = — -, \>p(5) (15) 

<Pd(\) - o 

Z {5 \x) = 1 + 6 [ W is) (y)dy, (16) 



where we recall that p(5) is solution to the Lundberg equation (pr)(\) = 6. Then from Expression (4) p. 19 of ffity 
one has that 

E[e- ST »] = Z^{b) - -L-WW(b). (17) 

Just to be clear on notations, we emphasize that [[20]] deals with spectrally negative processes. To apply it here 
(hence to obtain Expressions (|l5|), ( |l6| ) and (|l7|)), we thus need to consider hitting time of of process D t :— 
—D t starting from Dq — b. In particular, Laplace exponent ip(.) of D t as defined in Expression (2) of J2"(i| | by 
E[e A£)f ] = does coincide with function <£>£>(.), and $(<5) = sup{A > 0| i/;(X) = <5}, also defined in [Q, 

coincides with p(8). 



Remark 2.2 (scale function regularity) A necessary condition for function W^ s ' defined in the Proposition 2.3 



to be differentiate is that {D t , t > 0} has unbounded variation, which is the case here since it has a Gaussian 
component (i.e. a > 0). In fact it is shown in /^/ the stronger fact that a > implies that is twice 

differentiable. 

Remark 2.3 (boundary value of scale function) Still in the present case where process {Dt, t > 0} has un- 
bounded variation, we have that (0) = by Lemma 8.6 p.222 of[^. 



As a complement to (|17|), it is interesting to note that Remark 3 of [|20J gives an explicit expression of the joint 
Laplace transform of (XJ,, Dx b )- 

The approach in Proposition 2.3 has however a cost, which is that a Laplace Transform inversion of (|15|) is 



required to obtain the scale function. However recent results have been found concerning expression of ' in 
particular cases, see Hubalek and Kyprianou [|l8| as well as Egami and Yamazaki [|l4| in the case where {Gt, t> 
0} is a compound Poisson process with jumps foll owi ng p hase-type distribution. In fact the following result 
combines both approaches given in Propositions 22 and 23, and theoretically gives a closed form expression of 
scale function W^ s ' of any spectrally positive Levy process: 

Proposition 2.4 Scale function uniquely defined by Laplace transform satisfies the following first order 
differential equation 

W^'(x) - P (6)W^(x) = -?(p-f^g*i°(6,.)*h'(8,.)(5,x) := H(S,x) (18) 

fe=0 

where g(S, .) is given by f [/2[ ) and h'(S, .) is derivative ofh(S, .) given in with w = 1, i.e. 

h'(S,y) = -[-2fi/a 2 + p{8)]e-^l° 2+p ^v + 4 / e^- 2 ^ 2 +^ x -y^Q{x)dx 



- \-2^/a 2 + p(8)}— / e -[-^/° 2 +p(S)](ys) / e -*W*-')Q( x )dxd8. (19) 
° Jo Js 

Thus (x) has the following explicit expression 

W {5) {x)= f e- p{5){x - y) H{8,y)dy. (20) 



Proof: Differential equation (|l8]) simply come s fro m ( JTt| ) that one differentiates with respect to b (which is possible 
since is differentiable in light of Remark |2~2| ), using expression E[e~ l5Tb ] = <fr w (S, b) where penalty function 
w(.) is identically equal to 1, and finally using Expression (|l4|). Note that differentiation of ( |l4| ) is done by using 
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features derivative of function h(6, .). 



Since by Remark 2.3 one has that W^(0) = 0, Equation ( p0| ) is then obtained by solving the standard first order 
differential equation (|l8]). □ 

Note however that Formula (|20|) requires to compute the infinite series appearing in (fl8|), which in practice may 
not be handy. However, since such scale functions are important in the theory of Levy processes (in particular, 

and [| for determining quantities related to first passage times of 



these functions will be useful in Sections 2.3 



reflected processes and last passage times), any expression can be considered as welcome. 

Remark 2.4 Asymptotic behaviour ofTb as b — > +oo may be obtained through Roynette etal /p7|/. More precisely, 
it was proved that (T& + b/ ip' D (0)) / y/b converges in distribution to an Af(Q, — <p^, (0)/ (0) 3 ) distribution. One 

can also compute from /^/ asymptotic behaviour of triplet (t^ + b/ip' D (0))/Vb, Dx b — b, b — DT b -^jthat we 
did not include here but that involve technical expressions. 

2.2 Examples 

We illustrate the previous study with examples and review some famous examples related to degradation models. 

Pure gamma process Here we assume that o = and that {Gt , t > 0} is a gamma process with shape parameter 
a and scale parameter £. We recall that its Levy exponent and Levy measure are given by 

ip G {u) = <Pd(u) = -alog(l + u/£) 
vjj(dx) — Q(dx) — x~ 1 e~z adx. 

Considering this special case into the generalized Lundberg equation, it follows that this equation has no positive 
solution. It appears that the presence of the perturbation in the degradation model is important for applying the 



result obtained by Tsai and Wilmott [ |29| ] as we did in Proposition 2.2. However, in this first special case, the 
degradation process reduces to a pure stationary gamma process and so {Dt, t > 0} has increasing paths. It 
follows that: 

Vt>0, P[T b > t]=P[D t <b}. 

Consequently it is sufficient to study the distribution of Dt for any t > 0. The hitting time distribution was already 
given for instance p.517 of Park and Padgett 



Proposition 2.5 (Park and Padgett Q23|]) The cumulative distribution function (cdf) ofTb is: 

T(cd,b/0 



Vf > , F(t) 



T(at) 



where •) is the upper incomplete Gamma function. The probability distribution function (pdf) ofTb is, for any 
t > 0: 

fit) = a (*(«*) - log (|) ) 2^ + J -^ ) Q) ^ 2 F 2[ «t, at; at + 1, at + 1; - 6/0, 

where ^ is the di- gamma function (or logarithmic derivative of the Gamma function), •) = T(-) — r(-, •) is the 
lower incomplete Gamma function and 2 F 2 the generalized hypergeometric function of order (2, 2). 

It has been proved (see [jj]] or Section 5 of for instance) that Tf, has an increasing failure rate. 

Perturbed gamma Process Statistical inference in a perturbed gamma process has been studied in Q8|l u sing only 
degradation data. However sometimes both degradation data and failure time data are available (see J22| ] for such 
problem for a related model). In addition, from parameters estimation (based on degradation data for instance), 
one can obtain an estimation of the failure time distribution. Hence the distribution of TJ, can be of interest. In 
that case, {G t , t > 0} is a gamma process with shape parameter a and scale parameter £. We recall that Levy 
exponent and Levy measure of process {D t , t > 0} are then given by 



- a log(l + u/0 + \u 2 a 2 
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Thus, Proposition 2.2 gives joint distribution of (Tb, Dt,,-, Dt,,) through expression of <fi(S,b) where w(x) := 
f™w(x,y - x)^p-dy and g(5,y) = £ $ e - p ^ y -^ e'''™ 1 ^' ^^dxds, w(.,.) being an arbitrary 
bounded function. Also note that from Remark one has thanks to [|7| the Central Limit Theorem 



Finally, expression of the scale function is then given by ( |2C| ) with ipr>(-) and Q(.) defined in (j2ll). This will come 
in handy in Section || 

Brownian motion with positive drift We consider the case where Gt = lit, i.e. {D t , t > 0} is a Brownian 
motion with drift. In such case, the distribution of the hitting time of the constant boundary b is known and is 
called the inverse Gaussian distribution. Its pdf is given by: 



Vi>0,/(<) = ^=exp 



(b - ptf 



Proof of this result is generally based on the symmetric principle full-filled by the Brownian motion when \i = 0, 
or can be showed with martingale methods in the case p, > 0. Alternatively the pdf can be obtained by inverting 
the Laplace transform of Tb'. 

„-5T h i ( (IS ~ p)b 



m = E[e-^] = exp a z J , (22) 

with 7,5 = \/ p? + 28a 2 . Note that the expression of this Laplace transform is standard and can be found e.g. 
in Expression (38) p. 212 of [ |l2| ] (see also [y, page 19). Also note that ( ^2| ) is compatible with Expression 
dl4j). Indeed in the context of Proposition [^2 we have here 5 = and h = 0, thus ( |l4l ) reduces to <f>(6, b) = 
e -[-2fi/<r +p(S)]y Yyjjgj-g satisfies (^) •<=>■ = ^-p(S) 2 — pp(S) — S, giving the exact same expression (^2|). 
Expression of the scale function for this case is then given e.g. p. 121 in [n8(| by 



W {5) {x) = ; 2 e- px ' a2 sinh (^-^/2Sa 2 +p 2 ) = Ag-W* 2 sinh (JU^ . (23) 
y/25a 2 + A< 2 ® ' is Vcr 2 J 



Note that there seems to be a small mistake in JTs] | of expression of (x) (where there are some /i's instead 
of /i 2 's), that we corrected here. As proved by Chhikara and Folks [|Io|, the failure rate of an inverse Gaussian 
distribution is non-monotone, but it is initially increasing and then decreasing. 

Perturbed compound Poisson process with phase-type jumps Let us suppose that {G t , t > 0} is a compound 
Poisson process of intensity A whose jumps are phase-type distributed with representation (to, a, T). Let t := 
— Tl where 1 is a column vector of which entries are equal to l's of appropriate dimension (see e.g. Chapter VIII 
of Asmussen [||] for an extensive account on such distributions). In that case cpr) is given by 

(p D (u) = ^u 2 a 2 + X(a(ul - T) _1 t - 1). 

Egami and Yamazaki [|l4j] give the expression of the Laplace transform E(e~ l5Tb ) by determining a closed formula 
for the scale functions W 5 and using results in Proposition 2.3. More precisely following [53], let us denote for all 
5 > the complex solutions (£i ; s)i (resp. (r)i)i) of Equation <Pd(u) = S (resp. 5/ (5 — (fn(u)) = 0), u € C. We 
suppose that the £j,<s's are distinct roots. We set 

Is := {*| M-6,a) = * and > 0}, 

Js := {i\ 5/(5-<p d (-th)) = 0andH(»fc) > 0}, 



^5 ( U ) 



On page 4 of Jl4| ] it is stated that Card(Z,5) = Csrd(J's) + 1 (this results in fact comes from Lemma 1 (1) of | 
so that <pj(oo) exists and is equal to 0. We then define 

(Ai,s)iez, s.t. ipj (u) - ipj (00) = tpj (u) = Ai,s j^T > 



QS '■= Ai,s£,i,5 
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Then Proposition 2.1 of [ |l4| | gives expression of the Laplace transform <p{8) — E(e _l5Tb ) = 
J2i£i s Ai t se~^ i - S ^ x ~ a ^a.nd Proposition 3.1 of Jl4| ] yields the following interesting and useful expression of the 
scale function 



2 a Ct,<5 



e p{S)x _ e -£i,6X 



(24) 



Furthermore, as pointed out in [[l4|], expressions are more complicated but available when roots £ 4 ^'s have 

multiplicity mi > 1. 

2. 3 Reflected processes 

The previous model may not be too realistic if we consider the Brownian motion as a means of modelling small 
repairs, as the degradation process {D t , t > 0} may then be negative. An alternative for this is to consider the 
reflected version of {D t , t > 0} defined in the following way 

Vt>0, D*:=D t - inf (D s A0). 

0<s<t 

The hitting time distribution T h * of {D^ , t > 0} jointly to the overshoot and undershoot pdf is given by the 
following theorem 

Theorem 2.1 Let us suppose that {Dt, t > 0} is non monotone, i.e. that a > 0. Let be defined by 
where we recall that p = p{5) is solution to the Lundberg equation <Pd(z) = S. Then 

E[e- 5T >'l {D * edy , D . 6<to} ] = v D {dz - y)f ( b S) (b,y)dy (25) 

1 b J b 

Proof: We apply results from Doney Jl3| ] and we write, following notations therein, X t := —D t , so that Levy 
measure of {X t , t > 0} is H(dx) := V£,{— <ix)and process Y(t) := svp Q<s<t (X s V 0) — X t is equal to D* t . 
Following terminology of [|3J, is the (5-scale function of {X t , t > 0} and is defined by (|l5]) with ip-D instead 
of tpjj. Remark 4 p. 14 of [ pj[ | gives expression (|5|) where f £ is given by Pistorius [ ^4| (see also Expression (15) 
in Theorem 1 of Jl3|]) with x := and a := 6, noting that function is differentiable by Remark Q| . □ 



Note again that Theorem 2.1 is especially interesting when function admits closed form expressions, as 
m 00]. For example in the case of a perturbed compound Poisson process with phase-type distributed jumps 



(and using the same notations as in Section 2.2) we have V£>{dz) = Aae Tz t and given by (Q) (of which 



derivative is easily available), which, plugged in (|25|), easily yields the Laplace transform of the corresponding 
hitting time T b * jointly to the overshoot and undershoot distribution. 

We now state a famous lemma that links distribution of D* t to the cumulative distribution function of T b for all 

b > 0: 

Lemma 2.1 We have for all bandt> 0, P(D'* > b) = P(T b < t). 

Proof: This is a simple consequence from e.g. Lemma 3.5 p. 74 of Kyprianou Jl^jthat implies that P(Dj > b) = 
P (sup 0<s<t D s > b) which in turn is equal to P(T b < t). □ 

3 Last-passage time as failure time 

We let Li, and L* h be the last passage times of processes {D t , t > 0} and {D't, t > 0} below level b defined as 
L b := sup{0 <u\D u < b} and L* b := sup{0 <u\D* u < b} 
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3. 1 General case 

Let us introduce the following bivariate measures U and Won [0, +oo) 2 through their double Laplace transforms 



e- as - f3x U(ds, dx) = P[ } ~. P . , Vf3 > p(a), / / e'^-^Uids, dx) = — -, VjS, a > 0. 



o Jo at-ip D {P) Jo J p{at) + P 

(26) 

Expressions ( p^ ) may be found in Expressions (12) and (13) of [0], or p. 154 and p. 170 in Chapter 6 of jl9| ] (note 
that the latter reference considers spectrally negative processes, hence roles for U and U are swapped therein). 
Furthermore, from (26) of ^ one has that Us(dx) := f™ Q e~ Ss U{ds, dx) has the expression 

U s (dx) = e~ p ^ x dx, (27) 

hence 14$([0, +oo)) = l/p(5). In the same spirit, we define U$(dx) := J^ e~ Ss U(ds,dx). (|26|) then reads that 
I^Lo e ~^ x M$(dx) = s-!pd(P) ^ or a ^ @ > ^ e tnen nave f°ll° wm g identity, that will be of interest later 



Lemma 3.1 One has 

U s {dx) = [~p(S)W (s \x) + W {5) \x)]dx. (28) 



Proof: From ( 15 ) we get the following 

lO rOO rOD 

e- f}x U s {dx)^-p{8) e~ fix W {5) {x)dx + f3e- px W {5) (x)dx 
=o Jo Jo 



(29) 



where (3 > p{6). We recall from Remark 2.3 that W^(0) = 0. As to behaviour at +oo of the scale function, 
we have, thanks to Lemma 8.4 p.222 of [|9j, relation W {s) (x) = e ^W^ Vn( - c)) (x), for any c e R such that 
5 — tpr)(c) > 0, where w( s ~v d ( c '> is a scale function defined under a different probability measure. By picking 
c = p{5) then one gets 5 — ipo (c) =0 and 

W {5) (x) = eM'Wfil) (x) (30) 

(see e.g. Second Remark p. 32 of [ j25j ] for this identity as well as details on this other probability measure). At 
the end of Proof of Corollary 8.9 p.227 of JI^], it is shown that Wjf^ (+°°) — — ^~Jo+) wnere <PD,p(S) (?) : — 
<PD{q + p(S)) - pd{p{5)) = ifiD{q + p(S)) - S, hence 

Thus in view of (0) = 0, © and ©, and since f3 > p(S), the following integration by parts makes sense: 

1 r too r°° r°° 

0e- p "WW(x)dx= -e- p *wW(x) + / e~ Px W^ {x)dx = + / e"^ (x)dx, (32) 



remembering that is indeed differentiable by Remark 2.2. Comparing Laplace transforms (E9|) and (32), we 



then obtain (28h. □ 



Let us also note that, according to Definition 6.4 p. 142 of [[190, the fact that a > entails that is regular for 
sets (— oo, 0) and (0, +oo) (in particular, Theorem 6.5 p. 142 of jl9| ] applies here). With that in mind, and since 
{Dt, t > 0} is spectrally positive and drifts to +oo, we may recall the following important recent result from 
Kyprianou et al 

Theorem 3.1 (Corollary 2 of Kyprianou, Pardo and Rivero [|2l]]) Let us define 

Doo := inf D u , D s = inf D„ := sup{s > 0| D s - D s = 0}, 

u>0 t<s 

D* := inf D„. V * := inf-f.s > il D„ — D* = 01. 
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Then distribution of(G^, Dqq, ^L b — Lb, Lb,D f L b — b,b — D^ b ^, DL b ~b) is given by the following identity for 
t,b,v > 0, s > r > 0, < y < b + v, w > u > 0: 

PfG^ E dr, -Doo € dv, 2? Lb - L b E dt, L b E ds, T^ Lb -bEdu,b- D Lb _ E dy, D Lb -bE dw] 

= U s ([0, +oo)y 1 U(dr, dv)U{ds -r,b + v- dy)U(dt, w - du)Q{dw + y), (33) 

where U S {[Q, +00))^ = p(0)from @. 

It is clear that distribution of (Lb, b — D^ b ^, Di Jb — b) may be theoretically obtained from this theorem. In fact, our 



goal is to propose expressions of this distribution that only involves quantities that were determined in Section 2. 1 
e.g. scale functions, which we saw can be available in many situations, as opposed to measures U and U appearing 
in (|3j) which, as seen in (p6|), are available only through double Laplace transforms. More precisely, we have the 
following results. 

Theorem 3.2 We have for all t > and a E M, 



¥{L b <t) = / E[Di]W(o-6)/u t (o)da (34) 

Jb 

\L b >t, D t E da) = [1 - E[£>i]W(a - b)]f Dt (a)da (35) 



where fa, (■) is density ofr.v. D t and W(.) — W^ S ' (.) defined in (J73|) with 8 — 0. Besides, for all 5 > 0, and for 
b > y > 0, w > 0, the Laplace transform of Lb jointly to density of the under and overshoot is given by 



E[e b \b-D Lb -ed y , D Lb -bedw}} — 



dy.[l - e' p ^ w ]Q{dw + y). (36) 



Let us compare results given in Theorem |3.2| with existing ones in the literature concerning last passage times 
of Levy processes. References [ p"l| | and [|5j give distributions of respectively last exit times and last exit times 
before an exponentially distributed time, in terms of their Laplace transform, for a similar class of Levy processes; 



however Theorem 3.2 is more adapted here as it directly gives its cdf jointly to the density of the overshoot, thus 
avoiding an inverse Laplace transform. As said before, the slight advantage of Formula ( |36j ) over ([33]) is that it 
only involves the scale function. 

Proof: Let us start by showing ( |34| ) and (|3~5|). Let t > 0. By definition of Lb we note that for all a > b 
event [Lb < t,D t E da] is equal to [D t E da, {D s } will not hit level b anymore after t]. Hence using the Markov 
property: 

P[L 6 <t,U t €do] = P a -6 [T = +oo] P [D t E da] 

where P a _f, [T — +oo] is the probability that process {D t , t > 0} starting from a — b will never hit and is given 
e.g. through Formula (4) p. 19 of @ by P a _ b [T = +oo] = E[£) 1 ]W(a - b) and P [D t E da] = f D Aa)da where 



/p. is the density of r.v. D t and W(.) = W^ \.) in (|T5|). By integrating a from b to +oo one gets (34). Equation 
(H|) stems from the basic equality ¥(L b >t,D t E da) = F(D t E da) - P(L b < t, D t E da). 
We now turn to (pq), and use Theorem |3.l| to this end. Since by Fubini theorem we have 



/■OO 

\e- SLb I {b -D Lb _edy, D Lb -bedw}} = / e- Ss P[L b E ds, b - D Lb _ E dy, D Lb -bE dw], 



's=0 

and in view of (f33h, one just needs to compute the following integral: 



poo poo p pin 

/ / / / e-'TlGn E dr, -D^ E dv, V Lb - L b E dt, L b E ds, 

J v =oJt=OJs>r>oJu=0 

D^L b - b E du,b - DL b - E dy, Dh b - b E dw] 

poo p pw poo 

= p(0) / ti(dr,dv)e- Ss U{ds-r,b + v-dy). / U(dt,w - du).Q(dw + y), (37) 
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which we strive to do now. The first integral in the righthandside of ( |3~7| ) verifies, 

U{dr, dv)e~ Ss U(ds -r,b + v-dy) 

rOO 

U(dr,dv) / e- 5s U(ds-r,b + v-dy) 

J s—r 

U(dr, dv)e- Sr U§(b + v-dy) 
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v=0J s>r>0 
OO />00 



v — (V r— 

oo poo 



v — CV r— 

OO /"GO 



e~ 5r U(dr, dv) W {sy (b + v-y)- p(5)W {S) {b + v-y) 
U 6 (dv) \w {sy (b + v-y) - p(S)W (s) (b + v - y)} dy 



v—QJ r— 

oo 



dy by Lemma 3.1 



v=0 

oo 



e- p(s)v dv \w {sy (b + v-y)- p(5)W {S) (b + v-y) 



dy by Mi). 



(38) 



Relation @ yields that e ~P^>W^ (b - y + v) = eP^ b -^wfl(b - y + v) which, from (0), tends to 



aP(S)(b-y 



as v — > +00. This justifies the following integration by parts: 



■u=0 



e -p(s)v W (sy r b + v _ y ) dv = y e - P (s)v W {8) t h + w _ y j ^ ^ + / /9 ( ( 5) e -p(' 5 )^ (<5) (6 + « - 

1 



u=0 



p(<5)(6-y) 



^(p(*)) 

which, inserted in ([38]), yields the following simplification 

U(dr, dv)e~ Ss U(ds -r,b + v-dy) 



W (s) (b-y)+ I p(S)e- p ^ v W {s) (b + v-y)dv, (39) 

Jv=0 



I v—OJ s>r>0 

The second integral in the righthandside of (p7|) verifies 



,p($)(b-v) 



W^(b-y) 



dy. (40) 



w poo 



U (dt, w — du) 



u=0 Jt=0 



Uo(w — du) 

) 

e - P ( )( W - U ) du hy m 



u=0 

1 , 



(41) 

□ 



Plugging @ and (Q into @ yields @. 

3.2 Examples 

We consider here some examples from those studied previously and for which last-passage time is relevant. 

Brownian motion with positive drift In the case where G t = pt, fi > and D t = Gt + aB t = pt + aB t , we 
have 



E[Di] = th 



W(a-b) = W {0 \a-b) = -e-^ a - b ^ a \mh(^-^-p) from 

p ' ' 

fn t (a) = 



,-(a- fit) 2 1 '(2<r 2 t) 



which, plugged in (|^) and (^5|), yields expression of the cdf t M> P[L& < t] as well as its cdf jointly to density of 
D t . Note that by deriving this expression of the cdf one obtains after some calculation the following density for Lb 

F[L b e dt] = -J^e-—^-dt, 



/2%t 

which agrees with the already known density of the last passage time of a Brownian motion with drift, see e.g. 
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Perturbed gamma process In the case where {Gt, t > 0} is a gamma process with shape parameter a and 
scale parameter £, densities of G t and aB t are given by fc t (u) = "f^jy ^jJt- and f a B t (u) = a ^ 2 -Kt e ~ U ^ ' 
We also recall that function H (5, x) defined in Proposition 2.4 has expression given in (|l8]) with characteristics of 
the gamma perturbed process being given by (plj). Hence a bit of calculation yields 



W(a - b) 



i>a—b 



-p(S)(a-b-y) 



H(6,y)dy, 







fD t (a) = fa t * fcrB t (a) 



+ 00 



aV2TrtT(at)£ at J 

a -a 2 /(2a 2 t) r+°° i f 9 /9„2 t \ \ 



o-V27rtr(at)^ Qt 
e 



u at-l e -«/5 e -(a-n) 2 /(2 CT 2 t) dUj 
(fit 







a 2 /(2er 2 t) _ z'+oo 

— == (aVi) 01 *- 1 / at x **- 1 e~ 



1J 1 / 2g^t 

2 X 2tr V ^ 



^ £ 2a ) x dx i x ■— u/(aVi), 



2-KT{at)t c 



D— at 



<j\ft a 



where T{s) = / °° e~~ t t s ~ 1 dt, s > 0, is the gamma function and D p (z) = r(- P ) e~ zx ~ x l 2 x~ p ~ 1 dx, p < 0, 
is the parabolic cylinder function (see (9.241 .2) p. 1064 of Jl7|]). These expressions, plugged in ( |3~5| ) and ([36|), yield 
expression of the cdf of Lb jointly to density of D t as well as the Laplace transform of Lb jointly to density of the 
over and undershoot. 

Perturbed compound Poisson process with phase-type distributed jumps In the case where {G t , t > 0} is a 

compound Poisson process with phase-type distributed jumps of parameters as in Section 2.2, we have, using same 
notations as in that section that density of shocks is equal to p(x) = ae xT t (see Theorem 1.5(b) p. 218 of P])and 



E[Di] = -aT" 1 !, 
W(a-b) = -^-J2 A *,° 



£i,0 



a 2 g a 



ieio 



fo t (a) 



p(0) + 6,o 

' r).i 



e p(0)(a-b) _ e _£ i>Q (o-5) 



from (E4j) with 6 = 0, 



-u I {ia t) xjjggg expressions, plugged in (BBband ([36]), yield expression of the cdf of 



where f aBt (u) = ^ 

Lb jointly to density of D t as well as the Laplace transform of Lb jointly to density of the over and undershoot. 

3. 3 Reflected processes 

As for the previous section dealing with first-passage time, we consider the last-passage time for the reflected 
version of perturbed increasing Levy process. 



Theorem 3.3 The Laplace transform of L\ is given by 



E 


' e SLl 


= E[D 1 ] f 






Jb 



where we recall that i 



, a) = E[e 



-8T a 



W'{a-b)<j>(5, a)da 
4> w (5, a) with w = 1. 



Proof: We start similarly as in the proof of Theorem 3.2 and let T an independent r.v. follow- 
ing an £{8). Event [L\ < T, D* T e da] is equal to [Z?J € da, {D* s } will not hit level b anymore after T). 
Since reflected process {D$, t > 0} behaves like the non reflected process {D t , t > 0} on event 
[{D*} will not hit level b anymore after T] for t > T, we have, for all a > b, and using the Markov property, 



TO T T * ^ T n* r- A„~l — TO . PTL — I TO \ 7~1* t~ A„~\ 
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where P a -b [To = +00] is the probability that process {D t , t > 0} starting from a — b will never hit and has 
expression E[Di]T / K(a — b), as observed in Proof of Theorem 3.2. Since W(z) = on z < 0, we have by Fubini 
theorem (and since W(.) is a differentiable function by Remark |2l2"| ), 



E 




-I 






J a 



• [LI <T,D^e da] = E[Di] / W{a - b)P [D* T e da] 



a—b 



a—b 



= E[Di]E[W(D?p - b)} 



E[Di]E 



a—b 



W\a-b)l {D , T>a} da 



From Lemma 2.1, we have that P[D^ > a] 
distribution. This yields the result. 



/■OO 

= E[£>i] / W'(a - b)P[D* T > a]da. 

J a—b 

V[T a < T] which is equal to <f>(6,a), as T follows an £{5) 



□ 



Again we emphasize that </>(<5, a) 



-ST a 



is available in practice either through series (14) in Proposition 



2.2, or through ( |17| ) in Proposition 2.3. Also note that proof of Theorem 3.3 implicitly yields the following side 
result. 



Proposition 3.1 Let T be an independent £(5) distributed r.v. Then for all a > b we have 

d 

P[L* b >T, fl^e da] = -[1 - E[Dx]W{a - b)]—4>(5,a)da. 



(43) 



Proof: As in showing @, we use the fact that P[L* b >T,D* T ^ da] = P[D* T e da] — P [LI <T,DJ,e da] as 
well as @ to derive that P[L* b > T, D* T £ da] = [1 - E[D 1 ]W{a - b)]P[D* T G da]. To obtain flwe just 
need to prove that r.v. admits a density given by P[D^ g da]/ da — — ^0(5, a). Indeed Lemma 2A yields 



that P[D* T > a] = P[T a < T] = E[e' ST "] = <j>(5,a), thus what remains to prove is that E[e- 5T "] = (j>{8,a) 



is differentiable with respect to a. This can be seen thanks to the convenient expression ( |17| ) that yields that 
differentiability property since function is a differentiable function by Remark 2.2 (and is obviously 
differentiable by @). □ 



4 A maintenance policy 

We now as an application consider the maintenance strategy described in Barker and Newby H]. Degradation of 
a certain component is modelled according to a process {X t , t > 0}. We suppose that, without maintenance, 
{Xt , t > 0} is a perturbed process with same parameters as {D t , t > 0}and that failure occurs at the last passage 
time L b of level b of the degradation process. 

Let us then consider the following maintenance rule. The component is inspected at times (f/j)i=i,2,... such that 
inter inspection time verifies Ui+i — Ui = m(Xu i+ ), where m(.) is some non increasing function. Let d : K — > K 
be some "maintenance function". On inspection at time Ui, one of the following actions is undertaken: 

• either the system did not fail in interval (Ui-\, Ui], in which case preventive maintenance occurs and degra- 
dation process evolves like {D t , t > 0} with initial condition Dq = d(x) up until time Ui+\, where x is 
degradation state at instant Ui— ; in other words one has X\j i = d(Xjj i ^), 

• or the system failed in interval U] in which case it is repaired and degradation process starts anew, 
i.e. evolves like {D t , t > 0} with initial condition Do = 0. 

We will suppose in this section that function d(.) is differentiable from K to R and bijective. Note that these two 
assumptions are not too stringent and can be relaxed, in which case expressions of distributions computed in this 
section would only be more complicated. 

We then define r.v. / as the first inspection after which system is reset, i.e. 



T — 



inf In r~ T^T I 4""' 1 iln \'£* nii/iii rl in ( T T . 
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This means that T* := Uj is a regeneration time for the degradation process. Process {X t , t > 0} then behaves 
like independent copies of {D t , t > 0} in intervals (Z7j, J7j+i] with possibly different initial states. Figure [j] 
shows a sample path of {X t , t > 0}, with failure in interval (C/5, l^and thus starting anew at time Uq with 
Xjj 6 = 0. Note that process {X t , t > 0} thus constructed is cadlag and such that, given its state at any instant Uk, 
{Xt, t > Uk} is independent from {X t , t <G [0; Uk)}, i.e. from its history before Uk. This can be written as 



X t , t > Uk 



x s 



e[0,U k ] £ 



x t , t > Uk 



X. 




U, U 2 U 3 Vt (7 5 H b u 6 = U, = T' 

Fig. 1: Sample path of degradation process {X t , t > 0}, with failure in (Us, Ue]. 



We also introduce the idle time A* which is the unavailability period of time during which component is down 
until next scheduled inspection: 

A* :=T* -H b e [0,Ui-Ui-x] 

where Hi, is the failure time of the component and then necessarily lies in [£//-i,£/j]. We are interested in 
quantities involving (possibly joined) distributions of /, T*, A* as well as the state of the degradation process at 
inspection times. For this purpose we introduce the following quantities: 

• A[x, dy) := P[Lb > m(x), d(D m / x \) £ dy\ Do — x] the distribution of the degradation process on inspec- 
tion after maintenance jointly to the fact that there was no failure before inspection, given that degradation 
process starts at x, 

• C(y) := P[£fc < m(y)\ Do = y], the probability that failure occurred before next inspection, given that 
degradation process starts at y, 

• C r (y, z) := P[m(y) — Lb > z\ m(y) > Lb, Dq = y], z £ [0, m(y)}, the survival function of the idle time 
given that degradation process starts at y. 

These three quantities are easily obtained: 
Proposition 4.1 We have the following expressions 

A{x,dy) = [1 - npdWidr^y) -b + x)] ^-ifr)] ^ d V> 



C(y) = / nDi]W{a-b + y)f m(y) {a)da, 

Jb-y 

1 r°° 

C r (y,z) = — r-y / E[D 1 ]W{a-b + y)f m(y) ^ z {a)da. 
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Proof: We recall that we supposed that d(.) is a one to one differentiable function out of practicality. Expression 
for A(x, dy) simply comes from ( [35| ) with t = m(x) and a simple change of variable a — dr 1 {y)anA remarking 
that last hitting time of level b of process {D t , t > 0} with Dq = x is the same in distribution as that of level b — x 
of process {D t , t > 0} with Dq = 0. Expression for C{y) is obtained from ( p4| ) with t = m{y) and b := b — y 
because of process starting from y. Finally expression for C r (y, z) comes from the fact that 



C r {y,z) 



P[m(y) - L b > z\ D = y] F[m(y) - L b > z\ D Q = y] 



F[m(y) > L b , \D = y] C{y) 
and using ( [34] ) with T = m{y) — z and b := b — y to obtain expression of V[m(y) — L b > z\ Do = y]. 



□ 



We may now state main results of this section that concern quantities of interest introduced at the beginning of 
the section. 

Theorem 4.1 Distribution of I jointly to the state of the degradation process just after inspection and preventive 
maintenance is given by 

P[J = i, X Vl £ dyi,..., X Vi _ x € dyi-i] = A(0,dyi) x A{ yi ,dy 2 ) X ... x A(yi- 2 , x C{ Vl -i). (44) 

Distribution of the idle time jointly to I and the state of the degradation process just after inspection and preventive 
maintenance is given by 

P[A* > z, I = i, X Ul+ e dyi, X Ui _ 1 E dy^i] = A(0, dyi)x A(y ly dy 2 )x ...x%_ 2 , dj/<-i)xC r (j/i_i, z). 

(45) 

Proof: The first probability is obtained by writing it in the form P [n^L^-Efc n Fi\ where 



E k = 
Ft = 



no failure in (U k -x; U k ], d(X Uk ) e dy k 
failure in ({7j_i; Ui] 



Since evolution of process X t in t £ [Ui, Ui+i) given Xjj i is independent from X t , t £ [0, Ui), we may write that 
probability in the following form 



*[I = i, X Ut e d yi , X Uz _, e dyi-!] = [] F [E k \ x Uk _ t = y k _{\ x P[F,| X Vl _ t = W-i] 



fc=l 



and conclude by the fact that by the stationary increment property we have F[E k \ Xjj k l = y k -i] — A{y k _\, dy k ) 
and ¥[Fi\ Xu i l — yt-i] — C(y,_i) in order to obtain ([S|). (^5|) is derived by similar arguments. □ 



Note that Theorem 4.1 yields other interesting quantities. For example the expected time before reparation 
jointly to the number of inspections/maintenances is obtained thanks to (Q) by 



E [T* 



(j/i.-.w-OeR*- 1 



A(0,dyx) x A(y 1 ,dy 2 ) x ... x A(y^ 2 , X C(y i ^ 1 ). 



Remark 4.1 (Case of the reflected process) It is possible to adapt the previous setting to the reflected process 
{D% , t > 0} and constructed a reflected degradation process {X*, t > 0} with inspection and maintenance by 
considering exponentially distributed inter-inspection times Ui+i — Ui of which conditional distribution given Xjj i 
is £'(l/m(Xy i )), instead of deterministic times, where m(.) is the same function as in the non reflected caseand 
again featuring a maintenance function d{). Results from Theorem ^. ^ as well as equality ( f^i| ) would yield similar 
expressions for A(x, dy), C (y) for exponentially distributed horizonand an equivalent of Theorem \4.l\ for such an 
insnpction stratepv could he obtained 
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